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L  INTRODUCTION 


A.  GOAL 

The  two  most  \videly  used  models  for  area  search  are  Exhaustive  Search 
(ES)  and  Random  Search  (RS). 

If  searcher  speed  is  V,  detection  range  is  R,  and  area  size  is  A,  then  ES 
predicts  that  the  time  of  initial  detection,  T,  is  uniformly  distributed  between 
0  and  A/2RV.  This  model  is  appropriate  when  the  searcher  attempts  to  cover 
A  in  an  efficient,  systematic  fashion  and  the  target  is  stationary  and  uniformly 
distributed  over  A. 

Though  RS  is  traditionally  derived  with  a  stationary  target  and  a  random 
searcher,  it  also  covers  the  scenario  of  a  stationary  searcher  and  a  random 
target.  Therefore,  if  the  searcher  is  stationary  and  the  target  is  moving  over  A 
in  a  random  fashion  with  speed  U,  then  RS  predicts  that  T  will  have  an 
exponential  distribution  with  mean  A/2RU. 

The  goal  of  this  thesis  is  to  produce  and  test  a  new  analytical  model  which 
is  a  hybrid  of  ES  and  RS.  Both  the  searcher  and  target  will  be  allowed  to 
move.  The  searcher's  motion  will  be  systematic  (as  in  ES)  ,  and  the  target's 
motion  will  be  random  (as  in  RS).  Both  ES  and  RS  will  be  special  cases  of  the 
new  model,  called  Area  Motion  Search  (AMS). 

B.  METHODOLOGY 

Following  the  development  in  Washburn's  Search  and  Detection, 
[Ref.  l;pp.  2-1  -  2-4]  and  Koopman's  Search  and  Screening  [Ref.  2]  a  Poisson 


process  model  will  be  proposed  for  the  detection  process.  That  is,  an 
instantaneous  detection  rate  y(i)  \  "  be  introduced.  Then  the  cumulative 
distribution  function  (cdO  of  T  will  be 

f  t 

F(t)  =  1  -  exp  -Jy(u)du 

The  analytical  contribution  of  the  thesis  is  the  development  of  this 
detection  rate  model.  The  model  will  b<  tested  against  a  discrete-event 
simulation. 

Detection  rate  mathematics  and  the  development  of  the  AMS  detection 
rate  models  are  explained  in  Chapter  II,  Model  Development. 

Chapter  III,  Simulation  Development,  details  the  modification  of 
LT  Richard  L.  Darden's  simulation.  Hounds  and  Hares.  After  validation,  this 
s  mulation  was  used  to  examine  the  AMS  model.  Simulation  results  were 
compared  with  AMS  using  several  detection  rate  functions  in  an  effort  to  find 
the  best  model.  This  proce.ss  is  described  in  Chapter  IV,  Dat  Analysis. 
Finally,  Chapter  V,  Conclusions  and  Recommendations,  summarizes  and 
critiques  the  results  of  the  data  analysis  chapter  and  presents  areas  that  could 
u.se  further  research. 


II.  MODEL  DEVELOPMENT 


This  chapter  develops  the  AMS  model.  The  Assumptions  section 
examines  the  assumptions  necessary  to  determine  the  detection  rate  function 
for  AMS.  They  agree  with  both  the  ES  and  RS  model  assumptions.  The 
Detection  Rate  Background  section  explains  detection  rate  theory  and  the 
Area  Motion  Search  Detection  Rate  Function  section  applies  detection  rate 
theory  to  the  AMS  scenario  to  determine  a  basic  AMS  model.  Various  model 
variations  are  discussed  in  the  Embellishments  section. 


A.  ASSUMPTIONS 

The  AMS  stochastic  model  is  based  on  assumptions  that  are  consistent 
with  the  assumptions  used  in  the  ES  and  RS  models  [Ref.  l:pp.  1-2,  2-2  and 
Ref.  21: 

•  The  target  is  in  random  motion  on  .search  area  A,  with  speed  U.  (In 
ES  the  target  s  speed  U  equals  zero.) 

•  The  searcher  performs  a  systematic,  non-overlapping,  exhaustive 
search  at  constant  speed  V,  over  area  A. 

•  The  searcher  detects  the  target  with  probability  1.0  v/hen  distance 
from  the  searcher  to  the  target  is  within  detection  range  R,  and  does 
not  detect  the  target  otherwise.  (The  simplest  versions  of  classical  ES 
and  RS  are  actually  based  on  a  line  sensor  with  width  2R.) 

Detection  events  during  non-overlapping  time  periods  are 
probabilistically  independent. 

The  last  assumption  is  the  least  intuitive  and  the  most  questionable.  The 
degree  to  which  simulation  results  match  the  model  predictions  will  be  a 
measure  of  the  suitability  of  this  assumption. 
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B.  DETECnON  RATE  DEVELOPMENT 


A  detection  rale,  Y(t),  multiplied  by  a  very  small  time  period.  At,  is  the 
probability  of  initial  detection  in  that  small  time  period,  assuming  the  target 
has  not  yet  been  detected.  Therefore,  if  T  is  the  time  of  initial  detection  with 
pdf  f(t)  and  cdf  F(t),  then 

P(t<T<t+At)  f(t)At 
Y(t)At  -  pa>t)  -  1  -  F(t) ' 


or. 


Since 


dF(t) 


dt 

equation  is 


"  1  -  F(t)  • 

=  f(t),  it  is  easily  verified  that  a  solution  to  this  differential 


/  • 


lu 


F(t)  =  1  -  exp  l-J  Y(u)di 

I  0 

With  this  cdf  for  detection  time  T,  the  mean  time  of  detection  is 


E[T]  = 


exp 

-j  Y(u)du 

k  0  J 

dt. 


For  RS,  the  probability  of  target  detection  during  time  interval  At,  given 

no  earlier  detection,  is  the  area  searched  during  At,  2RVAt,  divided  by  the 

total  area,  A.  Applying  the  detection  rate  definition, 

P(t<T<t+At)  2RV 
"  P(T>t)At  ■  A 
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Therefore  RS  detection  time,  T,  is  an  exponential  random  variable  with  mean 

A  ,  ,  , 

equal  to  1  / y  =  •  iRf  f-  2] 

In  ES  the  target  is  stationary  so  the  maximum  time  spent  searching  is  the 

A 

time  needed  to  cover  the  entire  target  region  once,  or  2^^-  E)uc  to  the  fact 

that  the  target  is  at  a  uniformly  distributed  position  somewhere  on  the  target 
plane,  the  time  to  discover  the  target  is  a  uniform  random  variable  over 
[0,  and  the  mean  time  to  detection  is 

E[T|  = 


For  a  more  complete  discussion  sec  Washburn  [Ref.  Tpp.  1-1  - 1-4]. 

Even  though  detection  rates  are  generally  not  used  to  model  ES,  ES  can  be 
mathematically  manipulated  into  a  detection  rate  problem.  In  this  model, 
the  probability  that  the  target  has  not  yet  been  found  is  the  total  area.  A, 
minus  the  area  surveyed,  2RVt,  divided  by  the  total  area. 


.  (i>T) 


1  -  F(t) 


A  -  2RVt 
A 


consequently,  after  differentiation  and  algebraic  manipulation,  the  pdf  for  T  is 

,  2RV 
f(t)  =  -r-  - 


Therefore  from  the  definition  of  detection  rate, 

2RV 


7(0  = 


A  - 2RVt 


and  tlie  mean  detection  time  is 
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A 

4KV’ 


This  agrees  with  the  more  conventiortal  method  used  to  derive  ES. 

The  preceding  exanrples  suggest  that  the  motion  search  problem  could  be 
solved  as  a  detection  rate  problem.  The  trick  is  to  determine  the  detection 
rate  function,  7(1). 


C  AREA  MOTION  SEARCH  DETECTION  RATE  FUNCTION 

For  search  models  with  a  target  uniformly  distributed  over  a  fixed  search 
area,  the  detection  rate  can  be  interpreted  as 

_ rate  of  area  search  at  t _ 

^  ~  area  where  target  can  be  located  at  t ' 

In  AMS  the  area  behind  the  searcher  is  swept  free  of  targets.  This  area  is 
known  as  the  target  free  area  (TFA)  and  is  illustrated  by  Figure  1.  The  TFA 
has  zero  target  density,  while  the  rest  of  the  search  area  has  a  uniform  target 
density  greater  than  zero.  Therefore,  the  area  where  the  target  can  be  located 
at  time  t  is  the  total  search  area  minus  the  current  TFA.  This  means  the  AMS 
detection  rate  is 

rate  of  area  search  at  t 
~  search  area  -  TFA{t) 
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1.  Cookie-Cutter  Sensor  Target  Free  Area 

For  a  cookie-cutter  sensor  the  TFA  is  an  ice  cream  cone  shaped 
region,  as  in  Figure  2.  This  shape  was  developed  by  considering  a  target  that 
just  avoids  detection,  travelling  perpendicular  to  the  searcher  s  course.  Two 
such  targets  on  either  side  of  the  searcher  meet  at  center  of  the  searcher's 
path.  In  the  meantime,  the  searcher  has  advanced  some  distance.  These 
motions  combine  to  define  the  TFA  cone. 

AMS  starts  as  the  searcher  sweeps  into  the  search  area.  Fime  zero 
occurs  when  the  detection  disk  is  first  fully  in  the  search  area.  As  AMS 
progresses,  the  rest  of  the  TFA  is  pulled  into  the  search  area. 
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This  is  illustrated  in  Figure  3.  The  cone  portion  of  the  TFA  is 

attached  to  the  detection  disk  at  the  point  on  the  disk  where  the  slope  equals 

V 

U/ V,  point  a  in  Figure  2.  This  occurs  at  angle  6,  where  0  =  arctanf^  ). 


U 


Figure  3.  TFA  Cone  Development 
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The  cone  is  fully  pulled  into  the  search  area  in  the  time  it  takes  a 
target  to  travel  to  the  center  of  the  searcher's  path  from  point  a,  t  =  — —  , 

minus  the  time  the  searcher  requires  to  move  from  point  a  to  point  b, 
R(1  -  cos9) 

- Y - .  Therefore,  the  mechanics  of  TFA  development  can  be  explained 

by  two  equations.  One  equation  describing  the  sweeping  of  the  full  TFA  into 
,  ^  .  RsinG  R(1  '  cos6) 

the  search  area,  for  time  t  <  — y —  - y -  ,  and  another  equation  after 

the  entire  TFA  is  in  the  search  area. 

Table  1  illustrates  the  areas  of  several  geometric  shapes  used  to 
derive  the  TFA  equations. 

Using  the  geometries  derived  in  Table  1  the  TFA  equations  for  the 
ice  cream  cone  shape  are: 

(  a\D2  ,  d2  .  Q  R{l-cos9)^^  Rsin6  R(l-cos0) 

(jt-0)R  +  R  sin0cos9  -  UV[  t  + - - - )  :  t  <  “Tr - - - 

R(l-cos0)^ 


TFA(t)=^  +2VtRsin0ft4 


(7t~0)R^  +  R^sin9cos0  +  -^R^sinS 


:  i> 


RsinG  R(l-cos0) 
"T  "  V 


V 

where  0  =  arctanf  y  )■ 


These  equations  must  incorporate  one  more  contingency.  Physically, 
the  TFA  can  never  be  greater  than  the  search  area.  As  U  approaches  zero  (as 
in  ES)  the  TFA  grows  larger  and  the  TFA  equations  will  model  a  TFA  that  is 
larger  than  the  search  area.  But,  in  reality,  no  more  TFA  will  be  drawn  into 
the  search  area  after  the  searcher  has  traversed  the  entire  area  once.  The  time 
require  foi  a  single  searcher  traverse  is  • 
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,  A  Rsin0  R(1  -  cos6)  ,  ,  .  . 

Therefore,  if  <  — ^(j — ‘  - "y - equation  is 

truncated  at  time  t  =  and  the  TFA  remains  constant  at  that  value  for  the 

rest  of  the  search. 


TABLE  1.  TFA  GEOMETRIES  AND  AREAS 


Geometry 


Area 


(K  -  e)R2 


R^sinGcosG 


(V/U)R2sin2e 


2Vt(Rsin9)  -  UVt2 


2.  Line-sensor  Target  Free  Area 

In  ES  and  RS  the  TFA  is  nR^  when  the  searcher  is  motionless.  Both 
classical  ES  and  RS  models  ignore  this  area,  implicitly  assuming  line-sensors 
rather  than  circular  cookie-cutter  sensors.  The  ice  cream  cone  shape 
developed  for  the  AMS  model  is  based  on  a  cookie-cutter  sensor, 
consequently  it  will  not  ignore  the  tiR2  TFA  and  not  simplify  to  the  ES  and  RS 
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TFAs  in  the  extremum.  An  AMS  TFA  based  on  a  line-sensor  will  simplify  to 
the  actual  ES  and  RS  TFAs.  Therefore,  a  line-sensor  AMS  TFA  was 
developed  and  shaped  as  a  triangle  with  width  2R  and  length  RV/U.  For  a 
line-sensor. 


TFA(t)  =  i 


2RVt  -  VUl^ 
VR‘/U 


t<  RAJ 
t  >  RA’. 


These  equations  are  modified  to  reflect  the  possibility  that  the 
analytical  TFA  is  larger  than  the  search  area,  as  in  ES.  This  leads  to 

r 


2RVt  -  VUr 


■  r 

t  <  mir| 


A 

2RV 


TFA(l)  =  < 


miru 


a\t 


VR  . 

4R^V 


•  r  R  A  , 

.  t_mnY^,,2RvJ 


3.  AMS  Search  Rate  and  Dynamic  Enhancement 

For  many  search  models,  the  rate  of  area  search  equals  the  sweep 
rate.  This  is  incorrect  for  AMS  because  it  does  not  account  for  target  motion. 
Target  motion  increases  the  rate  of  area  search. 

A  model  must  be  developed  to  dynamically  enhance  the  searcher 
sweep  rate  with  the  target  motion.  U  and  V  are  velocity  vectors  of  known 
length  (V  and  U)  and  uniformly  random  direction.  (Bold  letters  refei-  to 
vectors,  plain  letters  refer  to  scalar  length.)  Therefore,  the  angle  between  U 
and  V,  G,  is  uniformly  distributed  on  [0,2k).  The  vectoi  addition  of  U  and  V, 
referred  to  as  V®U,  also  varies  with  9.  This  is  illustrated  in  Figure  4. 
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Figure  4.  Vector  Addition  of  V0U 

The  median  of  I V0U  I  is  VV-  +  and  occurs  when  U  and  V  are  at 
right  angles.  This  is  because  G  is  uniformly  distributed  over  (0,2:1)  and 

I  V®U  I  is  greater  than  Vv^  +  when  0  is  (-:t/2,n/2)  and  less  when  0  is 
(7c/2,3:t/2).  The  mean  of  I  V®U  I  is  given  by  Wa.shburn  as 

E[IV©UI]  =  2(V+U)E(K)/n, 


where 


^  "  U  +  V 


and 

E(K)  =  J  Vl  -  K^sin^O  d(>. 
0 


E(K)  is  a  complete  elliptic  integral  of  the  second  kind  and  cannot  be  explicitly 
evaluated,  but  is  tabulated.  Therefore  E[ !  V0U  I  ]  can  also  be  evaluated  from 
tables,  and  ranges  from  1.0  to  1.27  times  the  maximum  of  U  or  V.  [Ref  l:pp.  6- 
1  -  6-4] 
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For  purposes  of  model  introduction,  both  dynamic  enhancement 
models  will  notated  by  V©U. 

4.  AMS  Line-sensor  Detection  Rate  and  Expected  Detection  Time 

Using  the  definition  from  the  beginning  of  the  chapter,  the  AMS 
line-sensor  detection  rate  is 


(V©U)2R 

:  t  <  min 

■  R 

A  ■ 

A 

-(2RVt-VUt^) 

.U’ 

2RV. 

< 

(V©U)2R 

(V®U)8rV 

■  R 

A  ■ 

max 

_  A  -  rH’AJ  ’ 

A^U  J 

:  t  >  min 

.  U 

2RV_ 

Notice  that  if  U  or  V  equals  zero  these  equations  reduce  to  the  ES 
and  RS  detection  rates,  respectively. 

I'he  cookie-cutter  sensor  detection  rate  is  determined  analogously. 
These  detection  rates  can  be  olueeed  into  the  detection  rate  cdf  to  determine  a 

«  Vi/Sr^ 

closed  form  solution.  The  mean  time  to  detection  for  the  line-sensor  model 
^  R  A 

with  jj  <  2"^  can  be  determined  by  integrating  1  -  cdf  from  zero  to  infinity. 


EiTi  = 

Jo  [  Jo  A- 


(V©U)(2R) 


2RVt  +  Vui^ 


dt 


f  fR/^'  (> 

expl  -  - 

'  Jo  A- 


(V©U)(2R) 


A  -  2RVt  +  VL’E 


dt 


dt 


-I 


I  (V©U)(2R)  ,  ^ 
"v—  dt 

AJ  A  r,2  V 


A-r- 


dt. 


Applying  the  integral  formula 
dx  2 


1 


-tan 


.  1  2aX+b 


ax  +bx+(.  ^  4ac-b^  V  4ac-: 


(valid  for  o-c  4ac) 


where 


1.3 


a  =  VU 
b  =  --2RV 
c  -  A 


yields  the  open  form  solution  for  expected  detection  time. 


where. 


Bot  =  2V VUA  --  (RV)^ 
Z  =  (VeU)2R. 


D.  EMBELLISHMENTS 


The  proposed  A.MS  detection  ra 


tcs  do  net  cover  dll  possible  sceiidno^. 


Tliey  can  be  embellished  to  reflect  various  possible  occurrences. 

A  possibility  exists  that  there  is  an  area  larger  than  the  TFA,  behind  the 
searcher,  with  a  target  density  less  than  that  of  the  rest  of  the  search  area.  This 
area  will  be  called  an  effective  target  free  area  (ETFA).  A  TFA  and  a  larger 
ETFA  are  illustrated  in  Figure  5. 

Both  edge  effects  of  the  search  area  and  the  fact  that  the  target  rarely 
travels  directly  perpendicular  to  the  searcher  can  cause  an  ETFA  to  develop 
that  is  larger  than  the  TFA.  This  leads  to  adding  a  constant  of  multiplication, 
Ru,  to  modify  the  TFA  to  cTeate  the  longer  ETFA.  For  example,  if  Ru  equals 
two,  the  ETFA  is  twice  the  TFA.  The  resulting  detection  rate  is 
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f  (V©U)2R 

•  t  ^  TV11  n 

A 

1  A 

< 

max 

-(2RVi-VUi^./Rj 

(V©U)2R  (V®U)8R,R^V‘ 

.  1  ^  II 111 L 

1 

;  t  >  mii^ 

.'Tr’2^. 

A 

A-R„R^VAj'  ’  A^U 

U  ’2RV. 

Another  problem  arises  because  the  AMS  detection  rates  do  not  account 
for  the  searcher  turning  at  corners.  As  the  searcher  turns,  it  sweeps  out  less 
area  due  to  the  rotation  involved  in  turning.  This  means  the  sweep  rate 
decreases.  One  way  to  account  for  this  effect  is  to  v'se  an  effective  velocity  for 
the  searcher,  V'eff,  instead  of  the  true  velocity,  V.  V^ff  will  be  less  than  V  and 
will  depend  on  the  search  pattern  and  search  area  geometry.  Vpff  is 
determined  by  dividing  the  area  covered  in  one  cycle  of  search  by  the  time  it 
takes  to  perform  a  search  cycle  and  by  the  searcher's  sweep  width. 
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Therefore,  various  permutations  of  the  detection  rate  model  must  be 
investigated: 

•  Median  versus  mean  dynamic  enhancement,  discussed  in  the 
previous  section. 

•  V  versus  Vgff  searcher  velocities. 

•  Larger  effective  target  free  area  versus  the  actual  geometrically 
possible  target  free  area. 

•  Line-sensor  versus  cookie-cutter  sensor. 

The  resulting  detection  rate  functions  will  be  compared  with  the 
Hounds  and  Hares  simulation  results  in  the  next  chapter.  Data  Analysis,  in 
an  attempt  to  conclusively  support  one  line-sensor  detection  rate  model  and 
check  this  model  against  the  equivalent  cookie-cutter  model. 
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ni.  SIMULATION  DEVELOPMENT  AND  VALIDATION 


LT  Richard  L.  Darden  developed  Hounds  and  Hares  in  the  fall  of  1991. 
The  simulation  was  written  in  MODSIM  11^’'^^,  a  modular,  object-oriented 
simulating  language.  After  adding  a  patrolling  searcher  and  modifying 
the  randomly  moving  target  the  model  was  empirically  validated 
based  on  the  RS  and  ES  models.  This  chapter  examines  important 
aspects  of  the  simulation’s  design.  It  explores  search  area  design, 
searcher  patrol  parameters,  the  mechanics  of  the  random  target  and 
validation  of  the  simulation. 

A.  DESIGN  CONSIDERATIONS 
1.  The  Random  Target 

A  truly  random  target  has  a  position  at  one  instant  that  is  completely 
independent  of  the  position  at  the  next  instant  in  time.  It  has  the  exact  same 
probability  of  existing  at  one  point  as  any  other  point.  A  stationary  searcher 
would  detect  the  target  in  the  same  average  time  regardless  of  the  searcher's 
position  on  the  fixed  search  area. 

This  condition  is  difficult  to  meet  in  a  simulation  and  does  not 
reiiect  reality.  A  target,  in  reality,  has  a  continuous  path.  Therefore  its 
position  is  not  independent  from  one  instant  to  the  next  instant.  The 
random  target  in  Hounds  and  Hares  follow's  a  continuous  path.  The  target  is 
generated  at  the  beginning  of  a  simulation  run  at  a  uniformly  random 
position.  It  commences  travel  cn  a  uniformly  random  course.  When  the 
target  hits  a  search  area  boundary  it  is  reflected  by  the  diffuse  reflection 


process  suggested  by  McNish  in  Effects  of  Uniform  Target  Density  on  Random 
Targets  [Ref.  3:p.  43]. 

Diffuse  reflection  is  based  on  Lambert's  Law.  Lambert's  Law  states 
that  light  reflects  from  a  perfectly  diffuse  reflector  according  to  IqCOsG,  where  lo 
is  the  incident  light  and  6,  (-ji:/2,7c/2),  is  measured  from  the  normal  of  the 
diffuse  reflecting  surface.  [Ref.  4] 

When  a  target  hits  an  area  boundary,  the  reflection  angle  is 
determined  by  selecting  a  standard  uniform  random  variable,  Z,  and  solving 

a  =  arcsin(2Z  -  1), 

where  a,  the  reflection  angle,  is  measured  from  the  normal  of  the  boundary. 

Specular  and  uniform  reflection  were  also  investigated. 

The  effect,  on  target  randomness,  of  changing  the  target's  course  at 
random  times,  called  course  change  times,  was  also  explored.  Course  change 
times  selected  from  exponential  distributions  were  investigated. 

Target  randomness  was  investigated  by  placing  a  stationary  searcher 
on  a  target  plane  containing  a  randomly  moving  target.  The  RS  model 
predicted  expected  detection  times  equal  to  100  (A=50,  U-10,  V=0,  R=1.25). 
Table  2  illustrates  various  reflection  characteristics  and  resulting  average 
detection  time  confidence  intervals. 

Table  2  shows  that  both  specular  reflection  with  course  change  times 
and  diffuse  reflection  with  no  course  change  times  yield  excellent  random 
targets.  The  diffuse  reflection  method  requires  no  predictive  model  to 
determine  the  best  course  change  time,  unlike  the  specular  method. 
Therefore,  in  the  simulation,  the  random  target  was  programmed  to  follow 
diffuse  reflection.  Figure  6  shows  the  detection  time  cdf  of  both  cases  from 
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the  data  used  for  Table  2  superimposed  on  the  ideal  cdf.  It  shows  how  close 
targets  using  both  reflection  methods  are  to  ideal  random  targets. 


TABLE  2.  TARGET  CHARACTERISTIC  COMPARISOiXS 

95%  Detection  Time 
Target  Characteristics  Confidence  Interval 


Uniform  reflection- 
no  random  course  changes 

113-128 

Specular  reflection- 
no  random  course  changes 

167-190 

Specular  reflection- 
mean  course  change  time  =  5 

96-109 

Specular  reflection 

mean  course  change  time  =  20 

87-98 

Diffuse  reflection- 
no  random  course  changes 

91-104 

Diffuse  reflection- 

mean  course  change  time  =  20 

104-118 

Diffu.se  reflection- 

mean  course  change  time  =  5 

96-109 

Target  density  in  the  corners  of  the  search 

target  density  in  the  center  of  the  search  area. 

area  was  slightly  less  than 

This  is  demonstrated  by 

comparing  the  mean  detection  time  for  a  stationary  detector  placed  at 
different  positions  on  the  target  plane.  The  RS  model  predicts  an  average 
detection  time  of  100  for  the  scenario  illustrated  on  the  contour  plot  of 
Figure  7.  The  contours  show  longer  average  detection  times  in  the  corners  of 
the  search  area  indicating  lower  target  density  in  the  corners. 
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CUMULATIVE  DETECTION  PROBABILITY 
0.2  0.4  0.6  0,8 


o 


TIME 


Figure  6.  Cumulative  Detection  Probability  vs  Time  for  an  RS  Scer\ario  with 

Different  Target  Characteristics 
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Figure  7.  Average  Detection  Time  Contour  Plot  for  an  RS  Scenario  with 

Predicted  Detection  Time  =  100 
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2.  The  Patrolling  Searcher  and  the  Search  Area 

The  search  area  was  designed  as  a  square  and  a  ladder  search  was 
used  as  the  searcher's  patrol  pattern.  A  ladder  search  consists  of  parallel  paths 
offset  by  track  spacing  equal  to  twice  the  detection  range  [Ref.  5:p.  55].  Both  the 
search  area  and  the  search  pattern  are  typical  of  real  world  conditions. 

Non-uniform  target  density,  shown  in  Figure  7,  can  distort  AMS 
results  if  the  searcher  spends  an  inordinate  amount  of  time  in  low  target 
density  areas  during  many  simulation  runs.  To  overcome  this  problem,  the 
searcher  starts  patrolling  from  random  positions  at  the  beginning  of  each 
simulation  run.  Consequently,  the  average  target  density  seen  by  the  searcher 
approaches  a  uniform  density  over  the  course  of  many  runs. 

B.  SIMULATION  VALIDATION 

Tlie  ES  and  RS  models  were  used  to  validate  the  simulation. 

Figures  illustrates  the  detection  times  of  1000  simulation  runs  of  RS, 
with  predicted  mean  detection  time  1300,  plotted  on  a  quantile-quantile  plot 
against  an  exponential  distribution  with  mean  1300.  The  results  fall  along  the 
quantile- quantile  plot's  x  =  y  line.  This  validates  the  RS  extreme  of  the 
simulation. 

Figure  9  illustrates  the  results  of  1000  simulation  runs  of  ES  plotted  on  a 
quantile-quantile  against  an  ideal  uniform  distribution  between  0  and  1332 
time  units.  The  range  of  this  distribution  was  based  on  the  effective  velocity 
of  the  searcher,  discussed  in  the  En'ibellishments  Section  of  the  previous 
chapter.  It  was  necessary  to  use  V^ff  because  V  itself  predicted  an  upper  limit 
of  1300  tim.e  units  for  the  conditions  of  this  model.  Thus  it  did  not  reflect  the 
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actual  results.  It  is  also  intuitively  correct  to  use  V^ff  in  this  case  because  V 
does  not  account  for  the  turns  the  searcher  makes  in  the  ladder  search.  This 
quantile-quantile  plot  validates  the  ES  extreme  of  the  simulation. 


0  400  800  1200 

TIME  OF  DETECTION 


Figure  9.  Quantile-Quantile  Plot  of  ActuaJ  \-5  Ideal  Detection  Times  for  an 

ES  Scenario 
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IV.  DATA  ANALYSIS 


Data  analysis  was  split  into  exploratory  analysis  and  advanced  analysis. 
Data  analysis  concentrated  primarily  on  the  line-sensor  AMS  model.  During 
exploratory  analysis  30  AMS  scenarios  were  simulated  to  determine  average 
detection  times  for  each  scenario.  The  simulation  results  were  compared 
with  detection  tim.e  predictions  from  various  line-sensor  models  AMS 
models.  The  initial  findings  were  thoroughly  explored  in  advanced  analysis. 

The  follov/ing  flowchart  illustrates  data  analysis: 

Exploratory  Data  Analysis 

Exploratory  Analysis  concludes  AMS  predicts 
detection  time  well  if  the  target's  density  is 

I  inifnrm  (tciWa  tho  TP  A 

. . 

Advanced  Data^Analysis 

Graphically  investigate  the  relationship  between 
TFA  length  and  search  area  length  while 
observing  the  target  density  outside  the  TFA. 

Pick  the  AMS  model  that  correlates  most  closely 
with  the  simulation  while  target  density  remains 
uniform  outside  the  TFA. 

Investigate  AMS  when  the  target  density  is  not 
uniform  outside  the  TFA. 

Observe  the  effect  of  a  cookie-cutter  sensor  vice 
a  line-sensor. 
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A.  EXPLORATORY  DATA  ANALYSIS 


During  exploratory  data  analysis  30  different  AMS  scenarios  were 
simulated  to  find  a  pattern  between  line-sensor  model  predictions  and 
simulated  detection  times.  The  30  cases  were  chosen  based  on  a  wide  range  of 
parameters.  Search  areas  were  squares  with  side  lengths  of  25,  50  and  100 
units,  target  speeds  were  1,  3,  5  and  10  units  per  time  interval  and  searcher 
speeds  were  5, 10,  and  20  units  per  time  intervals.  Detection  ranges  were  1,  3, 
5  and  10  units. 

Each  scenario  was  simulated  to  determine  a  95%  confidence  interval  for 
the  mean  detection  time.  The  width  of  the  confidence  interval  was  less  than 


20  percent  of  the  mean  detection  time. 

Predicted  detection  times  were  evaluated  by  the  basic  line-sensor 
detection  rate  model  using  botlt  meihods  of  dynamic  enhancement.  I  he  two 
detection  rates  are 


(VeU)2R 

A-(2RVt-VU:^) 


(VeU)2R  (V®U)8R^^ 
A  -  R^VAJ  ’  A^U 


.  f  R  AT 

•  r  R  A  1 


median  dynamic  enhancement  model:  V0U  =  v  -t- 

mean  dynamic  enhancement  model:  V©U  =  E[  1 V  ®  U  I  ]. 

The  predicted  times  from  both  models  were  compared  with  the  center  of  the 
detection  time  confidence  intervals  in  an  effort  to  find  a  relationship. 
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A  pattern  emerged.  Predicted  and  simulated  detection  times  were  similar 
if  the  time  required  for  the  searcher  to  traverse  a  length  of  the  search  area  was 
at  least  five  times  greater  than  the  time  necessary  for  a  target  to  travel  to  the 
center  of  the  searcher's  path  from  the  edge  of  the  searcher's  path.  Another 
way  to  look  at  this  ratio  is  the  length  of  the  .search  area,  L,  (L  is  the  square  root 
of  A  for  a  square  search  area)  compared  to  the  length  of  the  target  free  cone, 
RV/U.  In  symbolic  terms,  if 


then  predicted  and  simulated  detection  times  matched  well. 

It  should  be  noted  that  if  U  equals  zero  this  ratio  equals  zero  and  AMS 
simplifies  to  ES.  Predicted  and  simulated  detection  times  also  match  well  in 
this  situation  Data  analysis  approached  AMS  primarily  from  the  RS  extreme 
towards  the  ES  extreme. 

This  conclusion  makes  intuitive  sense.  The  greater  the  ratio  is,  the 

the  smaller  the  TEA  is  in  comparison  with  the  search  area.  Therefore,  it  is 

less  likely  that  the  TEA  interacts  with  search  area  edges.  Th,'  embellishments 

to  the  basic  model,  ETEA  and  effective  searcher  velocity,  are  largely  a  result  of 

edge  effects.  Consequently,  embellishments  to  the  model  are  unimportant  if 
,  .  LU 

the  ratio  is  greater  than  five. 


B.  ADVANCED  DATA  ANALYSIS 
Exploratory  data  suggested  that  if, 

LU 

RV  - 
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then  either  dynamic  enhancement  AMS  line-sensor  model  provides  accarate 

LU 

forecasts  of  detection  time.  If  the  ratio  is  below  five  other  factors  affect  the 
results. 

1.  ETFA  versus  TFA 

LU 

The  accuracy  of  the  basic  model  when  g-y  is  greater  than  five  leads  to 

the  conclusion  that  the  ETFA  approaches  the  TFA  for  these  scenarios.  To 

LU 

investigate,  a  scenario  with  greater  than  five  was  selected.  The  scenario 

was  simulated  many  times  (~1000)  while  the  searcher  was  started  in  a  corner 

of  the  search  area.  Each  simulation  was  stopped  when  the  searcher  reached 

the  center  of  the  search  area.  The  position  of  undetected  targets  was  plotted. 

The  geometric  TFA  was  overlaid.  Figure  10  illustrates  for  a  scenario  with  a 

ratio  of  8.33.  The  geometric  TFA  is  drawn  and  the  target  density  is  roughly 

constant  outside  the  TFA.  The  ETFA  equals  the  TFA.  These  results  are 

isistent  with  results  from  other  scenarios  with  ratios  greater  than  five. 

The  results  of  Figure  10  lead  to  an  interesting  question.  Is  the  ETFA 

>  iFA  for  low  ratios?  Graphs  created  by  the  previous  method  indicate  yes. 

One  ■  the  more  striking  is  for  a  scenario  with  a  ratio  equal  to  1.67,  illustrated 

in  Fiqure  11.  The  density  of  targets  outside  the  geometric  TFA  in  this  graph  is 

not  uniform.  'Fhere  is  a  lower  density  swath  that  corresponds  to  an  ETFA 

LU 

longer  than  TFA.  This  effect  disappears  as  the  ratio  is  increased. 

2.  The  Most  Accurate  Line-Sensor  AMS  Model 

LU 


py  greater  than  tive  were 

LU 


200  simulation  scenarios  with 
investigated.  They  were  divided  into  ten  groups  of  20  scenarios  with  = 
10,20,. ..,100.  The  cases  were  determined  by  randomly  selecting  A  between  25 
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NAUTICAL  MILES 


and  100,  U  and  V  between  1  and  20,  setting  the  ratio  constant  and  solving 
for  R. 


Figure  10. 


Target  Density  and  TFA  for 


LU 

RV 


=  8.33 


1 
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Figure  11,  Target  Density  and  TFA  for 

■  r 
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The  simulation  results  were  compared  with  four  detection  times 
predicted  by  four  different  line-sensor  detection  rate  models  consisting  of 
permutations  of  the  two  dynamic  enhancement  models  and  the  searcher  true 
and  effective  velocities.  Symbolically,  V©U  in  the  standard  line-sensor 
detection  rate  formula  was  replaced  by 

E[IV©UIL  EdVeffeUI],  and  • 

The  results  were  compared  by  determining  the  relative  error 
between  each  predicted  and  simulated  detection  time.  The  relative  error  of 
each  prediction  was  averaged  into  an  average  relative  error  for  each 
predictive  model.  The  results  are  illustrated  by  the  box  plot  in  Figure  12. 

The  mean  dynamic  enhancement  (V©U  =  E[  I  V  ©  U  I  ])  detection  rate 
model  provided  the  least  average  error  at  the  2%  significance  level  according 
to  Fisher's  least  significant  difference  method.  [Ref.  6]  This  model's 
predictions  had  an  average  percent  error  between  3.8%  and  4.8%  at  the  957c 
confidence  level.  Figure  13  illustrates  the  results  of  this  model  plotted  against 
actual  detection  time.  Figure  14  plots  the  residuals  of  Figure  13. 

3.  Effective  Searcher  Speed 

The  effective  searcher  speed  embellishment,  V^ff,  is  designed  to 
account  for  a  searcher  sweeping  less  area  in  a  turn,  leading  momentarily  to  a 
reduced  detection  rate.  It  is  interesting  to  examine  why  the  actual  searcher 
speed,  V,  provides  a  better  model  of  the  simulation  than  Vpff,  as  shown  in 
Figure  12.  One  reason  may  oe  that  V  provides  a  better  detection  rate  for  the 
great  majority  of  search  time  so  a  V^ff  based  detection  rate  is  unnecessary. 
The  physical  Hounds  and  Hares  simulation  characteristics  might  also  cause 
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REL^nVE  ERROR 

0.05  0.10  0.15  0.20 


this  effect.  The  reduced  target  density  at  the  edges  of  the  simulation  search 
area  may  already  compensate  for  the  turning  searcher  making  a  Veff 
correction  redundant. 


E[V+U]  (V^+u^)^  E[Vjfj  +  U]  (Vj;rr^+U2)^ 

MODEL 


Figure  12.  Box  Plot  of  Relative  Error  between  Simulation  and  AMS  Model 

Predictions 
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Figure  13. 
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Figure  Error  between  Predicted  and  Simulated  Detection  Times  for 


Scenarios  with  >10 


V 


LU 

ETFA  versus  TFA  for  Small  Scenarios. 


The  basic  AMS  model  does  not  successfully  predict  detection  times 
LU 

for  scenarios  with  small  ratio:.  This  is  illustrated  by  Figures  15  and  16.  In 

Figure  15  simulated  detection  times  are  plotted  against  predicted  detection 

LU 

times  for  scenarios  with  ratios  between  one  and  five.  Figure  16  is  a  plot  of 

the  resulting  residuals.  The  poor  correlation  between  predicted  and 

simulated  detection  times  is  the  result  of  large  TFAs  interacting  with  search 

area  edges  to  create  non-uniform  target  densities  outside  the  TFAs. 

LU 

Scenarios  with  smaller  ratios  were  investigated  to  see  if  the 

ETFA  embellishment  helps  predict  accurate  detection  times.  50  random 

LU 

scenarios  were  developed  for  each  ratio  m  the  set  (1,2,...,!  0).  The  detection 

rate  model  with  the  ETFA  embellishment  was  applied  to  the  scenarios. 

Each  scenario  was  examined  with  ETFAs  that  were  multiples  of  the 

TFA.  A  sensitivity  range  was  determined  for  the  TFA  multiples  based  on 

average  error  between  predicted  and  simulated  detection  bmes.  The  valley  in 

Figure  17  shows  the  ranges  of  multiples,  Ru,  that  yield  the  least  average  error 
,  ,  LU  LU 

for  each  ^  scenario.  Figure  17  shows  that  as  the  ratio  gets  smaller,  down 

to  about  one,  the  optimal  Ru  becomes  greater.  This  is  due  to  the  increasing 

interaction  of  the  TFA  with  search  area  edges. 

LU 

Figure  17  shows  that  for  large  j^y  ratios,  the  TFA  is  relatively  small, 

and  the  model  is  less  sensitive  to  increasing  ETFA.  In  fact,  the  results  in  the 

LU 

Figure  17  show  that  ETFA  effects  can  be  ignored  for  ratios  as  low  as  j^-y  >  5. 
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Figure  16.  Error  betweeen  Predicted  and  Simulated  Detection  Times  lor 

LU 

■  Scenario.s  with  p-y  <  10 
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Figure  17.  3-D  Contour  Plot  of  Av  rage  Error  vs  ETFA  Multiples 

■  r 
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Figures  18  and  19  plot  the  same  scenarios  as  Figures  15  and  16,  while 
applying  an  ETFA  correction  to  the  AMS  model.  The  results  are  vastly 
improved. 

LU 

In  order  to  examine  AMS  when  one,  50  random  scenarios 
LU 

were  simulated  for  each  ratio  in  the  set  (0.2,0.3,0.5,0.7,0.8). 

The  results  of  scenarios  with  ratios  below  0.8  are  nearly 

indistinguishable  from  ES  predictions  (as  if  U  =  0).  These  scenarios  reflect  the 

ES  extreme  of  AMS  because  the  target  has  little  probability  of  crossing  the 

searcher's  long  tail.  In  fact,  an  AMS  prediction  with  an  ETFA  embellishment 

is  only  slightly  better  than  an  ES  prediction  in  this  region. 

An  ETFA  embellishment  was  applied  to  these  scenarios.  Figure  20 

interpolates  and  plots  the  results.  Figure  20  shows  that  the  optimal  R-j 
LU 

decreases  as  decreases  below  one  because  these  scenarios  approach 

exhaustive  search.  The  optimal  Ru  forces  the  ETFA  to  equal  the  TFA  for  the 
equivalent  ES  scenario  with  U  equal  to  zero.  The  optimal  Ru  and  the  TFA 
tend  to  follow  the  rough  relationship 

(TFAtruc)  -  TFAes  ■ 

As  TFAtrue  becomes  larger,  Ru  becomes  smaller  to  maintain  constant  TFAes- 
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DETECTION  TIME 


Figure  18.  Corrected  Predicted  vs  Simulated  Detection  Times  for 

LU 

Scenarios  with  j^<  10 
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Figure  20.  Ru  and  Error  vs.  LU/RV  ratio 

5.  Line-sensor  versus  Cookie-Cutter  Sensor 

The  results  of  the  cookie-cutter  sensor  model  were  indistinguishable 

LU 

from  the  line  sensor  model  for  ratios  as  low  as  >  5.  This  is  expected 


because  in  higher  ratio  scenarios  the  model  loses  sensitivity  to  TFA,  as 

illustrated  in  Figure  17.  Therefore  the  slightly  larger  shape  of  the  cookie- 

cutter  TFA  will  not  affect  the  results.  Differences  betv/een  the  models  in 
LU 

lower  regions  are  visible,  requiring  different  ETFA  corrections  for  the 


two  models. 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 


This  thesis  attempts  to  shed  light  on  an  interesting  area  of  search  theory, 
area  motion  search.  While  investigating  this  phenomenon,  many  subtopics 
were  touched  upon:  the  characteristics  of  a  random  target  and  the 
relationship  of  edge  effects  on  search  models  are  just  a  few.  Many 
conclusions  can  be  drawn  from  this  thesis,  but  all  must  be  treated  skeptically 
as  they  are  only  as  good  as  the  supporting  simulation.  In  short,  this  is  a 
dynamic  area  of  operations  research  that  could  benefit  from  much  more 
study. 

A.  CONCLUSIONS 

The  area  motion  search  problem  can  be  treated  as  a  detection  rate 
problem.  Line-sensor  and  cookie-cutter  sensor  detection  rate  models  both 
predict  area  motion  search  well  when  the  target  is  uniformly  distributed  over 
the  search  area.  The  line-sensor  area  motion  search  detection  rate  reduces  to 
the  classical  exhaustive  search  detection  rate  and  random  search  detection 
rate  when  the  target  and  searcher  velocities  are  zero,  respectively.  The  line¬ 
sensor  detection  rate  is 


(V©U)2R 

t  <  min 

A  1 

A 

-(2RVt-VUl^) 

L  U’2RVj 

Y(t)  =  ■< 

(V0U)2R 

(V©U)8R\' 

t  >  min 

r  R 

A  1 

max 

< 

.  A  -  Rnw  ’ 

A^U 

.  U’2RvJ 

A  =  search  area,  V  =  searcher  velocity,  U  =  target  velocity,  R  =  detection  range 


If  the  target  is  bounded  by  the  search  area,  edge  effects  can  prevent 
uniform  target  distribution  over  this  area.  Edge  effects  can  Le  ignored  for  a 
square  search  area  when  the  length  of  the  search  area  is  on  the  order  of  five 
times  greater  than  the  length  of  the  target  free  area,  or 
LU 

^  5,  L  =  length  of  a  search  a  .‘ea  edge. 


When  edge  effects  cannot  be  ignored,  the  line-sensor  detection  rate  can  be 
modified  by  expanding  the  taigel  free  area  to  equal  an  effective  target  free  area 
to  accurately  predict  average  detection  times.  The  effective  target  free  area 
accounts  for  reductions  in  target  density  in  certain  parts  of  the  search  area  due 
to  the  interaction  of  the  true  target  free  area  with  the  edges  of  the  search  area. 

To  expand  the  target  free  area,  the  factor  Ru  was  introduced  into  the 
existing  line-sensor  detection  rate  to  create  a  new  detection  rate. 
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For  example,  if  Ry  equals  two,  the  effective  target  free  area  is  twice  the 
geometric  target  free  area. 

LU 

Figure  20  illustrates  optimum  Ru  multiples  versus  for  a  square  search 

LU 

area  and  a  ladder  search  pattern.  Figure  20  shows  Ru  increasing  as 

LU 

increases  from  zero  to  one,  Ry  peaking  near  =  one  and  Ry  decreasing  as 
LU 

increases  beyond  one. 
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Ru  decreases  as  increases  above  one  because  the  target  free  area 

becomes  smaller  in  relation  to  the  search  area.  This  makes  edge  effects  less 
prominent,  and  sntaller  effective  target  free  areas  are  required, 

Ru  increases  as  increases  from  zero  to  one,  because  these  scenarios 

approach  exhaustive  search.  The  optimal  Ru  multiple  forces  the  effective 
target  free  area  to  roughly  equal  the  target  free  area  for  the  equivalent 
exhaustive  search  scenario  (target  velocity  equals  zero),  according  to 

(target  free  areatruc)  *  ~  target  free  areaexhaustive  search- 

Therefore,  as  target  free  areatrue  becomes  smaller,  Ru  becomes  larger  to 
maintain  constant  target  free  areaexhaustive  search- 

B.  RECOMMENDATIONS  FOR  FURTHER  RESEARCH 

•1^1  ••i*  t  f  f  » y 

iiuS  ineSis  suggests  many  areas  ror  runner  researen,  inciuaing; 

•  Determine  the  optimal  real  world  search  pattern  for  a  searcher  in 
AMS. 

•  Attempt  to  duplicate  the  results  on  different  search  area  geometries 
with  different  search  patterns. 

•  Investigate  the  best  random  target  characteristics  for  simulation. 

•  Investigate  edge  effects  more  closely. 

•  Investigate  the  affects  of  introducing  a  smart  target  to  the  model. 

•  Attempt  to  validate  the  model  in  sea  tests. 

•  Introduce  the  area  motion  search  into  low-fidelity  combat  models. 
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